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Abstract 

The objective of present research is to examine the thermal radiation effect in three-dimensional mixed convection flow of 
viscoelastic fluid. The boundary layer analysis has been discussed for flow by an exponentially stretching surface with 
convective conditions. The resulting partial differential equations are reduced into a system of nonlinear ordinary 
differential equations using appropriate transformations. The series solutions are developed through a modern technique 
known as the homotopy analysis method. The convergent expressions of velocity components and temperature are 
derived. The solutions obtained are dependent on seven sundry parameters including the viscoelastic parameter, mixed 
convection parameter, ratio parameter, temperature exponent, Prandtl number, Biot number and radiation parameter. A 
systematic study is performed to analyze the impacts of these influential parameters on the velocity and temperature, the 
skin friction coefficients and the local Nusselt number. It is observed that mixed convection parameter in momentum and 
thermal boundary layers has opposite role. Thermal boundary layer is found to decrease when ratio parameter, Prandtl 
number and temperature exponent are increased. Local Nusselt number is increasing function of viscoelastic parameter and 
Biot number. Radiation parameter on the Nusselt number has opposite effects when compared with viscoelastic parameter. 
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Introduction 

Analysis of non-Newtonian fluids is an active area of research 
for the last few years. Such fluids represent many industrially 
important fluids including certain oils, shampoos, paints, blood at 
low shear rate, cosmetic products, polymers, body fluids, colloidal 
fluids, suspension fluids, pasta, ice cream, ice, mud, dough floor 
etc. In many fields such as food industry, drilling operations and 
bioengineering, the fluids, either synthetic or natural, are mixtures 
of different stuffs such as water, particle, oils, red cells and other 
long chain molecules. Such combination imparts strong rheolog- 
ical properties to the resulting liquids. The dynamic viscosity in 
non-Newtonian materials varies non-linearly with the shear rate; 
elasticity is felt through elongational effects and time-dependent 
effects. The fluids in these situations have been treated as 
viscoelastic fluids. Further, all the non-Newtonian fluids in nature 
cannot be predicted by single constitutive equation. Hence all the 
contributors in the field are using different models of non- 
Newtonian fluids in their theoretical and experimental studies (see 
[1-11] and several refs. therein). The boundary layer flows of non- 
Newtonian fluids in the presence of heat transfer have special 
importance because of practical engineering applications such as 
food processing and oil recovery. Especially the stretching flows in 
this direction are prominent in polymer extrusion, glass fiber and 
paper production, plastic films, metal extrusion and many others. 



After the pioneering works of Sakiadis [12] and Crane [13], 
numerous works have been presented for two-dimensional 
boundary layer flow of viscous and non-Newtonian fluids over a 
surface subject to linear and power law stretching velocities (see 
some recent studies [14-21]). It has been noted by Gupta and 
Gupta [22] that stretching mechanism in all realistic situations is 
not linear. For instance the stretching is not linear in plastic and 
paper production industries. Besides these the flow and heat 
transfer by an exponentially stretching surface has been studied by 
Magyari and Keller [23]. In this attempt the two-dimensional flow 
of an incompressible viscous fluid is considered. The solutions of 
laminar boundary layer equations describing heat and flow in a 
quiescent fluid driven by an exponentially permeable stretching 
surface are numerically analyzed by Elbashbashy [24] . Al- Odat et 
al. [25] numerically discussed the thermal boundary layer on an 
exponentially stretching surface with an exponential temperature 
distribution. Here magnetohydrodynamic flow is addressed. 
Nadeem and Lee [26] presented the steady boundary layer flow 
of nanofluid over an exponential stretching surface. Sajid and 
Hayat [27] examined the thermal radiation effect in the boundary 
layer flow and heat transfer of a viscous fluid. The flow is caused 
by an exponentially stretching sheet. The thermal radiation effect 
in steady hydromagnetic mixed convection flow of viscous 
incompressible fluid past an exponentially stretching sheet is 
examined by El-Aziz and Nabil [28]. Pal [29] carried out an 
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analysis to describe mixed convection heat transfer in the 
boundary layer flow on an exponentially stretching continuous 
surface with an exponential temperature. Here analysis is given in 
the presence of magnetic field, viscous dissipation and internal 
heat generation/absorption. Khan and Sanajayand [30] investi- 
gated the heat and mass transfer effects of viscoelastic boundary 
layer flow over an exponentially stretching sheet in presence of 
viscous dissipation and chemical reaction. Bhattacharyya [31] 
numerically investigated the heat transfer boundary layer flow 
over an exponentially shrinking sheet. Shooting method is 
implemented here. Recently, Mukhopadhyay et al. [32] dealt 
with the boundary layer flow and heat transfer of a non- 
Newtonian fluid over an exponentially stretching permeable 
surface. Mustafa et al. [33] studied the boundary layer flow of 
nanofluid over an exponentially stretching sheet with convective 
boundary conditions. Flow and heat transfer for three-dimensional 
viscous flow over an exponentially stretching surface is discussed 
by Liu et al. [34]. Bhattacharyya et al. [35] studied the effects of 
thermal radiation in the flow of micropolar fluid past a porous 
shrinking sheet with heat transfer. The transient free convection 
interaction with thermal radiation of an absorbing emitting fluid 
along moving vertical permeable plate is discussed by Makinde 
[36]. Hayat et al. [37] considered a two-dimensional mixed 
convection boundary layer MHD stagnation point flow through a 
porous medium bounded by a stretching vertical plate with 
thermal radiation. 

Literature survey indicates that the published studies about 
three-dimensional flow by an exponentially stretching surface are 
still scarce. To our knowledge, there is only one recent study by 
Liu et al. [34] which describes the three-dimensional boundary 
layer flow of a viscous fluid over an exponentially stretching 
surface. Thus motivation of present research is to venture further 
in the regime of three-dimensional mixed convection flow of 
viscoelastic fluid over an exponentially stretching surface with 
thermal radiation. The surface possess the convective type heat 
condition. No doubt the thermal radiation effects are significant in 
many environmental and scientific developments, for instance, in 
aeronautics, fire research, heating and cooling of channels, etc. It 
is found that radiative transport is often comparable and hence 
associated with that of convective heat transfer in several real- 
world applications. Therefore it is of great worth to the researchers 
to study combined radiative and convective flow and heat transfer 
aspects. Moreover, the skin friction coefficients for three-dimen- 
sional viscoelastic fluid have been computed which has not yet 
been available in the literature. This paper is structured into the 
following fashion. Section two consists of mathematical formula- 
tion and definitions of physical quantities of interest. Convergent 
series solutions of the involved nonlinear systems are developed in 
section three. The solutions in this section are developed by 
homotopy analysis method (HAM) [38-45] . Section four comprises 
discussion with respect to seven pertinent parameters involved in 
the solutions of velocity components and temperature. Section five 
syntheses the main observations. 
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Figure 1. Geometry of Problem. 

doi:10.1371/joumal.pone.0090038.g001 

The governing boundary layer equations for steady three- 
dimensional flow of viscoelastic fluid can be put into the forms (see 
Nazar and Latip [11]): 
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Mathematical Modelling 

We consider three dimensional mixed convection boundary 
layer flow of second grade fluid passing an exponentially stretching 
surface. The surface coincides with the plane z = 0 and the flow is 
confined in the region z>0. The surface also possess the 
convective boundary condition. Influence of thermal radiation 
through Rosseland's approximation is taken into account. Flow 
configuration is given below in Fig. 1. 



where w, v and w are the velocity components in the x — , y — 
and z — directions respectively, is the material fluid parameter, 
/I is the dynamic viscosity, v = (/i/p) is the kinematic viscosity, T is 
the fluid temperature, p is the fluid density, g is the gravitational 
acceleration, is thermal expansion coefficient of temperature, 
Cp is the specific heat, k is the thermal conductivity and qr the 
radiative heat flux. Note that w-momentum equation vanishes by 
applying boundary layer assumptions (see Schlichting [46]). 



PLOS ONE I www.plosone.org 



2 



March 2014 | Volume 9 | Issue 3 | e90038 



Three-Dimensional Mixed Convection Flow 



f"(0) g"(0).^'(0) 



-A"r-0.5,Pr- 1.2, 

Jv ■ U.J, K ■ A ■ V.m 
1 1 1 I ■ 1 ■■ 1 1 1 1 1 1 1 1 




01 




-10 -Of -0« -04 

m 








o: 
















-10^ 






1 if 







Figure 2. ^ — curvesforthefunctions/(^),g(^)and0(^). 
doi:10.1371/journal.pone.0090038.g002 



By using the Rosseland approximation, the radiative heat flux qr 
is given by 



qr-- 



(5) 



The boundary conditions can be expressed as 



u=U^,v= V^,w = 0, -k^=h(Tf-T)Mz = 0, 
oz 



Where c^is the Stefan-Boltzmann constant and kg the mean 
absorption coefficient. By using the Rosseland approximation, the 
present analysis is limited to optically thick fluids. If the 
temperature differences are sufficiently small then Eq. (5) can be 
linearized by expanding into the Taylor series about T^^, 
which after neglecting higher order terms takes the ft)rm: 



(8) 



where subscript w corresponds to the wall condition, his the 
thermal conductivity, Tf is the hot fluid temperature, his the heat 
transfer coefficient and Too is the free stream temperature. 

The velocities and temperature are taken in the following forms: 



T'^=4TiT-3Tt. 



(6) 



x+y x+y A(x+y) 

U„ = Uoe—, V„ = Foe—, T„. = Tf=T^ + Tqc-^ (9) 



By using Ens. (5) and (6), Eq. (4) reduces to • u- u tz ^ ^ r • r i j 

J ^ ^ \ f \ ^ \ J m which (7o, '^o arc the constants, L is the reierence length and 

^is the temperature exponent. 

dT dT dT k d^T IGasT^^d^T The mathematical analysis of the problem is simplified by using 

^~dx^^~dy^^~d^ ^ 'pcpj? ~ 3kepCp dz^ ' ^""^ the transformations (Liu et al. [34]): 



Table 1. Convergence of series solutions for different order of approximations when ^ =0.1, A = 0.2, Pr =1.2, a = 0.2, 
hf=-0.5fig=-0.6 and he=-0.1. 





Order of aproximations 


1 


5 


10 


15 


20 


25 


-f (0) 


1.06111 


1.02482 


1.02609 


1.02623 


1.02618 


1.02618 


-g" (0) 


0.544444 


0.548057 


0.548092 


0.548043 


0.548053 


0.548053 


-e' (0) 


0.317778 


0.305581 


0.305729 


0.305744 


0.305738 


0.305738 



doi:10.1371/joumal.pone.0090038.t001 
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Figure 3. Influence of a. on the velocity f'(r]) 

doi:10.1371/journal.pone.0090038.g003 



Figure 5. Influence of X on the velocity f'(rj). 

doi:10.1371/journal.pone.0090038.g005 
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/ = 0, g = 0,/ =hg' = oifi'=-y{\- em at ^ = 0, (14) 



Incompressibility condition is now clearly satisfied whereas Eqs. 
(2)-(7) give 



f"' + if+g)f"-2if'+g')f' + 
k{ 6/"'/' + (3."-3/" + .."')/" ^^^^^^^ 

V + (V + ing")f"' -if+g+ ng')/'" . 



(11) 



f'^0,g' ^0,0^0 asjy^oo 



(15) 



in which^ is the viscoelastic parameter, a is the ratio parameter, 
Pris the Prandtl number, Gfj^is the local Grashof number, i?is the 
radiation parameter, ^is the temperature exponent, y is the Biot 
number, Re^^ is the local Reynold number, k is the mixed 
convection parameter and prime denotes the differentiation with 
respect to f]. These can be defined as 



g"'+(f+g)g"-^if'+g')g'+ 

6g"V + (3/"-3g" + >/"'k" 



K 



-(4f'+2nf")g"'-(f+g+nf')g" 



(12) 



a - 0.2, R - 0.3, y - A - 0.4, X -0.5, Pr - 0.7 




a = 0.2, A =0.5, A = 0.4, K ^ R = 0.3, Pr = 0.7 



y- 0.1, 0.4, 0.7,1.0 




Figure 4. Influence of K on the velocity /'(^). 

doi:10.1371/joumal.pone.0090038.g004 



Figure 6. Influence of y on the velocity f'(f]). 

doi:10.1371/joumal.pone.0090038.g006 
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Figure 7. Influence of R on the velocity f'(r]). 

doi:1 0.1 371 /journal. pone.0090038.g007 



Figure 8. Influence of oc on the velocity g'(r]). 

doi:10.1371/journal.pone.0090038.g008 
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by 



The skin-friction coefficients in the x and y directions are given 



Re 



-1/2 



Qy = (^) [g" + K(-{f+g)g'" + 5{f' + g')g" + 2fT + 2g'g'X^,. (20) 



Further the local Nusselt number has the form 



Nu^ \,r\J''' --j{^Ti^4R)m. (21) 
k{T,,-T^)lx L\2 J 3 



Cfy- 



(17) 



Series Solutions 

The initial guesses and auxiliary linear operators in the desired 
HAM solutions are 



where 



/o(;7) = (l-e-''),go('?) = o((l-e-''),0oW=exp(-^) (22) 



du - 
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(18) 



z = 0 



By using Eq. (18) in Eq. (17) the non-dimensional forms of skin 
friction coefficients are as follows: 



Re 



-1/2 



r+K(-{f+ g)r + 5{f' +g')f" + 2/7" + 2g'g")]., ( 19) 
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Figure 9. Influence of K on the velocity g\r]). 

doi:10.1371/journal.pone.0090038.g009 



PLOS ONE I www.plosone.org 



5 



March 2014 | Volume 9 | Issue 3 | e90038 



Three-Dimensional Mixed Convection Flow 





a - 0.2, R = 0.3, y = .\ = 0.4, A =0.5, Pr = 0.7 



K = 0.0, 0.5, 1.0, 1.5 



Figure 10. Influence of a. on the temperature 9(r]). 

doi:10.1371/journal.pone.0090038.g010 



Figure 11. Influence of K on the temperature 6(r]). 

doi:10.1371/journal.pone.0090038.g011 



subject to the properties 



Lf(Ci + Cse^ + €36-^) = 0, Lg(C4 + Cse"^ + Cee'"^) = 0, 
LeiC,e^^Cse-^) = 0 



(24) 



in which C/ (/ = 1 — 8) are the arbitrary constants, Lf,Lg and Lq 
are the linear operators and /o(^),go(^) and ^o(^) are the initial 
guesses. 

Following the idea in ref. [38] the zeroth order deformation 
problems are 

(1 -p)Lf [fir,; p) -Mn)] =/'^/N/ [fin; p), gin; p)] , (25) 



(\-p)Lg[gin; P)-gQin)] =Ph^g [fin; P\g(n; p)] , (26) 



mation equations deform from the initial guesses /o(^),^o(^) and 
^o(^)to the exact solutions/(/7), g(^)and 6(r])of the original 
nonlinear differential equations. Such kind of continuous variation 
is called deformation in topology and that is why the Eqs. (26-28) 
are called the zeroth order deformation equations. The values of 
the nonlinear operators are given below: 



+ K 



drj drf" 



5^2 



\ 



fin,p)+gin,p)+n- 



dr] dr]^ j drj^ 

Sg(ri,p)\ s'^fin^p) 



dr] 



(30) 



(1 -p)Le [)ir, p)-eoifi)] =pmo [fim pXm pi kfi,P)] , (27) 



/(O; p) = 0,f{0;p)=\,f(ca;p) = 0, ^(0; = 0, 

^■(0; p) = a, g'irx> ; = 0, 9(0, p) = - y [1 - 9(0, p)], (28) 

9(-x),p) = 0, 

For p = 0 and p=l one has 



fir, 0)=Mrihm O)=go(rih 

9(fi,0) = 9o(fi),andf(m 1) =/('/), 
gir,; l)=g(^), %,!) = %). 



(29) 



Note that when p increases from 0 to 1 then f(f],p), g(f],p) and 
6(rj,p) vary from/o(^), go(rj) and ^o(^) to f(rj), g(f/)and^(^). So 
as the embedding parameter pE[0, 1] increases from 0 to 1, the 
solutions f(rj] p),g(r]] p) and 0{r]] p) of the zeroth order defor- 



0.3,a-r-A-0.4,Pr-0.7 




Figure 12. Influence of X on the temperature 0{r]). 

doi:1 0.1 371/journal.pone.0090038.g01 2 
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K - R - 0.3, a - A - 0.4, X - 0.5, Pr - 0.7 



y = 0.1, 0.4, 0.7, 1.0 





K = R = 0.3, a = y = 0.4, X = 0.5, Pr = 0.7 



A = 0.0, 0.3, 0.6, 0.9 



Figure 13. Influence of y on the temperature 9(r]). 

doi:1 0.1 371 /journal.pone.0090038.g01 3 



Figure 14. Influence of A on the temperature 9(r]). 

doi:10.1371/journal.pone.0090038.g014 
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Sg{n,p) 
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+ K 
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dr] dt]^ j drj^ 
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drj J drf 



0(;7,;7) = (?oWX^0,„('7K',0,„W = 



1 5"-%; p) 



ml dp^ 



p = 0 



, (35) 



where the convergence of above series strongly depends upon 
X (31) hffig^nd ho. Considering that fif, fig and ^gare chosen in such a 
manner that Eqs. (33)-(35) converge at />= 1 then 



(36) 



g(n)=go(n)+ ^gmin). 



(37) 



Ng[e{n,P)J(n,p),g{n,p)] =(i + 



df] 



(32) 



(38) 



Here ^/,^gand ho are the non-zero auxiliary parameters and 
N/,Ngand the nonlinear operators. Taylor series expansion 



1 ^J{n;p) 



ml dp^ 



, (33) 



gin-.p)=gdn)+ ^gm(ri)p"\gmiti)-- 



1 d"'g(r,;p) 



ml dp^ 



(34) 



p^O 



X - 0.5, Pr - 0. 




Figure 15. Influence of R on the temperature 9(r]). 

doi:1 0.1 371/joumal.pone.0090038.g01 5 
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0.4 



K - R - 0.3, a - 0.2, A - y - 0.4, A - 0.5 




Figure16. Influence of Pr on the temperature 6(f]) 

doi:10.1371/journal.pone.0090038.g016 

The corresponding problems at mth order deformations satisfy 
Lf\fUfl)-xJ,n-i(fl)]=fif'Rl'(fl), (39) 

Lg[g„,{r,)-x,„g,„-M]=figR'^iri), (40) 

Lg[dUfl)-X,ne,n-M]=fteKM. (41) 



Rg(n)=g„- l(l)-2Y,Sm-t -kSk -^^gm-\-kfk 

k-0 k-0 

m— 1 

^{f„,-X-kgk+gm-\-kgk) 



k = 0 
( 



k=0 k=0 k=0 

m — l m—l m—l 

k=0 k=0 k=0 

m-l ^^^^ m-l ^^^^ m-l 

- E gm-l-kg'k - J2 fm-l-kgk - E nfm-\-kSk 
k=0 k=0 k=0 



(44) 



+ FTk=o"'-\0^_,_,fy+0:„_,_,g,) 

E (45) 

m—l 
k = 0 

The mth order deformation problems have the solutions 

f,„(n)=f,:il)+Ci+C2e'' + Qe-'', (46) 



/„,(0) =U0) =U^) = 0, g,„iO)=gJO)=gJ^) = 0, 

(42) 

„(0)-ye,„(0) = e„(a)) = 0. 



gUfl)=g:„(n) + C4 + Cse" + Qe-", (47) 

f^™('?) = f^;„('?)+C7e'' + C8e-'', (48) 
where the special solutions are fj*, and 6^. 

Convergence Analysis 

We recall that the series (36-38) contain the auxiliary 
parameters fif, fig and ho. These parameters are useful to adjust 
and control the convergence of homotopic solutions. Hence the 
h— curves are sketched at 15^^^ order of approximations in order 
to determine the suitable ranges for hf, hg and ho. Fig. 2 denotes 
that the range of admissible values of hf, hg and ho are 
-0.7</z/<-0.2, -0.7<^^<-0.1 and -0.8 < -0.2. 
Table 1 shows that the series solutions converge in the whole 
region of ^ when hf= — 0.5 fig = —0.6 and ho= —0.7. 

Table 2. Comparative values of —f"(0), —g'XO) and /(oo)+g(oc) for different values a when Ki=X = y = R = 0. 



KJ^in)=f:-M-^Y.f'rn-l-kfk-^Y.Sn. 



+ ^(fm-\-kf'k+gm-\-kfk) 
k = 0 



/ m-l m-l 

' ^Y.fm-1-kfk^^Y.gm-l-k 
/t = 0 yt = 0 



/t = 0 
m-l 



■ E ng,n-l-kfk^^ E gm-l-kfk +2 E ngm-l-kfk 
k=0 k=0 k=Q 



- Y.fm-1-kfk - E gm-l-kfk - E ngm-l-kfk 
k=0 k=Q k=0 



(43) 



-2xe, 







Liu et al. [34] 








Present results 






a 


-f " (0) 


-g" (0) 


f(co)+g(co) 


-f" (0) 


-g" (0) 


f(co)+g(co) 


0.0 


1.28180856 


0 


0.90564383 


1.28181 


0 


0.90564 


0.50 


1 .56988846 


0.78494423 


1.10918263 


1.56989 


0.78494 


1.10918 


1.00 


1.81275105 


1.81275105 


1.28077378 


1.81275 


1.81275 


1.28077 
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Table 3. Values of skin friction coefficients for different values 
of K and a when X = y = 0.5, R = 03, Pr = 1.2 and A = 0.2. 



K 


a 






0.0 


0.5 


4.95289 


4.37363 


0.2 




5.16586 


3.97055 


0.3 




5.42622 


3.96130 


0.3 


0.0 


3.72170 


1 .65409 




0.2 


4.30247 


2.34617 




0.5 


5.42622 


3.96130 
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Discussion of Results 

The effects of ratio parameter a, viscoelastic parameter K, 
mixed convection parameter A, Biot number y and radiation 
parameter R on the velocity component f'(r]) are shown in the 
Figs. 3-7. It is observed from Fig. 3 that velocity component f'(f]) 
and thermal boundary layer thickness are decreasing functions of 
ratio parameter a. This is due to the fact that with the increase of 
ratio parameter a, the x-component of velocity coefficient 
decreases which leads to a decrease in both the momentum 
boundary layer and velocity component /'(vj). Fig. 4 illustrates the 
influence of viscoelastic parameter K on the velocity component 
f'{y{). It is clear that both the boundary layer and velocity 
component f'{j]) increase when the viscoelastic parameter 
increases. Influence of mixed convection parameter k on the 
velocity component f'{rj) is analyzed in Fig. 5. Increase in mixed 
convection parameter A shows an increase in velocity component 
f'ij]). This is due to the fact that the buoyancy forces are much 
more effective rather than the viscous forces. Effects of Biot 
number y and the radiation parameter R on the velocity 
component f'(r]) can be predicted from Figs. 6 and 7. These 
Figs, depict that the influences of y and R on both the velocity 
component f\f]) and thermal boundary layer thickness are similar 
i.e. there is increase in these quantities. Figs. 8 and 9 illustrate the 
variations of ratio parameter a and viscoelastic parameter K on 
the velocity component g'(r]). Variation of ratio parameter ais 
analyzed in Fig. 8. Through comparative study with Fig. 3 it is 
noted that fXf]) decreases while g'(^) increases when a increases. 
Physically, when a increases from zero, the lateral surface starts 
moving in y-direction and thus the velocity component g'ifj) 
increases and the velocity component f'(r]) decreases. Fig. 9 is 
plotted to see the variation of viscoelastic parameter K on the 
velocity component g'ifj). It is found that both the velocity 
component g'{r]) and momentum boundary layer thicknesses are 
increasing functions of It is revealed from Figs. 4 and 9 that the 
effect of ^ on both the velocities are qualitatively similar. Figs. 10- 
1 6 are sketched to see the effects of ratio parameter a, viscoelastic 
parameter K, the temperature exponent A, Biot number y, mixed 
convection parameter A, Radiation parameter and Prandtl 
number Pr on the temperature 6(ri). Fig. 10 is drawn to see the 
impact of ratio parameter a on the temperature 6(r]). It is noted 
that the temperature 6(f]) and also the thermal boundary layer 
thickness decrease with increasinga. Variation of the viscoelastic 
parameter K on the temperature 6(r]) is shown in Fig. 1 1 . Here 
both the temperature and thermal boundary layer thickness are 
decreasing functions of K. Variation of mixed convection 
parameter X is analyzed in Fig. 12. It is seen that both the 
temperature 6(r]) and thermal boundary layer thickness are 



decreasing functions of mixed convection parameter A. Fig. 13 
presents the plots for the variation of Biot number y. Note that 
6(ri) increases when y increases. The thermal boundary layer 
thickness is also increasing function of y. It is also noted that the 
fluid temperature is zero when the Biot number vanishes. 
Influence of temperature exponent A is displayed in Fig. 14. It 
is found that both the temperature 6(f]) and thermal boundary 
layer thickness decrease when A is increased. Also both the 
temperature 6(r]) and thermal boundary layer thickness are 
increasing functions of thermal radiation parameter R (see 
Fig. 15). It is observed that an increase in R has the ability to 
increase the thermal boundary layer. It is due to the fact that when 
the thermal radiation parameter increases, the mean absorption 
coefficient kg will be decreased which in turn increases the 
divergence of the radiative heat flux. Hence the rate of radiative 
heat transfer to the fluid is increased and consequently the fluid 
temperature increases. Fig. 16 is plotted to see the effects of Pr on 
6{r]). It is noticed that both the temperature profile and thermal 
boundary layer thickness are decreasing functions of Pr. In fact 
when Pr increases then thermal diffusivity decreases. This 
indicates reduction in energy transfer ability and ultimate it results 
in the decrease of thermal boundary layer. 

Table 1 presents the numerical values of —f"(0), —g"(0) and 
— 6'(0) for different order of approximations when 
fif= — 0.5 fig = —0.6 and fio= —0.7. It is seen that the values of 
—f"(0) and —g"(0) converge from 20th order of deformations 
whereas the values of —0'(0) converge from 25th order 
approximations. Further, it is observed that we have to compute 
less deformations for the velocities in comparison to temperature 
for convergent series solutions. Table 2 includes the values for 



Table 4. Values of local Nusselt number —6'(0) for different 
values of the parameters K, a, A, R, A, Pr and y. 





K 1 


a y R 


Pr 


A 




0.0 0.5 


0.5 0.5 0.3 


1.2 


0.2 


0.297492 


0.3 








0.308234 


0.5 








0.311853 


0.2 0.0 








0.303062 


0.3 








0.304775 


0.5 








0.305738 


0.2 0.5 


0.0 






0.282007 




0.3 






0.297135 




0.5 






0.305738 




0.1 






0.0885730 




0.3 






0.216850 




0.5 






0.305738 


0.2 0.5 


0.5 0.5 0.0 






0.329701 




0.3 






0.305738 




0.5 






0.292750 


0.2 0.5 


0.5 0.5 0.3 


1.0 




0.292152 






1.2 




0.305738 






1.5 




0.321826 


0.2 0.5 


0.5 0.5 0.3 


1.2 


0.0 


0.288530 








0.2 


0.305738 








0.5 


0.325492 
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comparison of existing solutions with the previous available 
solutions in a limiting case when Ki = ^ = y = R = 0 and a varies. 
This Table presents an excellent agreement with the previous 
available solutions. Table 3 is computed to see the influences of 
viscoelastic parameter K and ratio parameter a on skin friction 
coefficients in the x and y directions. It is noted that K has quite 
opposite effect on skin friction coefficients while quite similar effect 
is seen within the increase of ratio parameter a. Table 4 examines 
the impact of viscoelastic parameter K, mixed convection 
parameter A, ratio parameter a, Biot number y, radiation 
parameter R, Prandtl number Pr and temperature exponent A 
on the local Nusselt number (rate of heat transfer at the wall). It is 
noted that the value of rate of heat transfer increases for larger 
viscoelastic parameter K, mixed convection parameter X, ratio 
parameter a, Biot number y, Prandtl number Pr and temperature 
exponent A while it decreases through an increase in radiation 
parameter R. 

Table 2. Comparative values of —f"(0), —g"(0) and 
/(oc) + g(oc) for different values a when Ki=l = y = R = 0. 

Conclusions 

Three-dimensional mixed convection flow of viscoelastic fluid 
over an exponentially stretching surface is analyzed in this study. 
The analysis is carried out in the presence of thermal radiation 
subject to convective boundary conditions. The main observations 
can be summarized as follows: 

• Influence of ratio parameter a on the velocities f'(f]) and g'(ri) 
is quite opposite. However the effect of viscoelastic parameter 
K on the velocities f'(fj) and g'ivj) is qualitatively similar. 

• Momentum boundary layer thickness increases for g'(f]) when 
ratio parameter a is large. Effect of a on f'(r]) is opposite to 
that of 

• Velocity component f'(f]) is increasing function of mixed 
convection parameter A. However 6(f]) decreases with an 
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